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ABSTRACT 

Aims. I deduced a three-dimensional sunspot model that is in full agreement with spectropolarimetric observations, in order to address 
the question of a possible penumbral heating process by the repetitive rise of hot flow channels. 

Methods. I performed inversions of spectropolarimetric data taken simultaneously in infrared (1.5yum) and visible (630 nm) spectral 
lines. I used two independent magnetic components inside each pixel to reproduce the irregular Stokes profiles in the penumbra. I 
studied the averaged and individual properties of the two components. By integrating the field inclination to the surface, I developed 
a three-dimensional model of the spot from inversion results without intrinsic height information. 

Results. I find that the Evershed flow is harbored by the weaker of the two field components. This component forms flow channels 
that show upstreams in the inner and mid penumbra, continue almost horizontally as slightly elevated loops throughout the penumbra, 
and finally bend down in the outer penumbra. I find several examples, where two or more flow channels are found along a radial cut 
from the umbra to the outer boundary of the spot. 

Conclusions. I find that a model of horizontal flow channels in a static background field is in good agreement with the observed 
spectra. The properties of the flow channels correspond very well to the Moving Tube simulations of Schlichenmaier et al. (1998). 
From the temporal evolution in intensity images and the properties of the flow channels in the inversion, I conclude that interchange 
convection of rising hot flux tubes in a thick penumbra still seems a possible mechanism for maintaining the penumbral energy 
balance. 
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1. Introduction 

Sunspots have been the first indicators of solar magnetic ac- 
tivity detected already centurie s ago in the Western World 
(Galilei et al. 1613; Galilei 1632), respectively, millennia ago in 
Asia dWittmann &Xulll987l). Their mag netic nature was only 
proven in the last centu r y by | Hald (1 19081) . Using spectroscopic 
observations, lEvershedl (1 19091) showed that sunspots exhibit a 
particular flow field in the penumbra, the brighter halo surround- 
ing the dark umbra: as soon as a sunspot was observed not di- 
rectly on disc center, spectral lines in the half of the spot facing 
disc center were seen to be blue shifted, whereas they were red 
shifted in the half facing the limb. This "Evershed flow" was ex- 
plained as the signature of a radial outflow, which had to be close 
to parallel to the solar surface. 

Despite its rather well defined spectroscopic properties, the 
Evershed flow however still eludes a generally accepted ex- 
planation. The reason is the fine-structure of the sunspot's 
penumbra that became visible when both the spatial reso- 
lution and the information content of observations was im- 
proved. Spectroscopic observations with high spatial resolu- 
tion show that the Evershed flow is organized in many small- 
scale radially oriented filaments, which carry the bulk of the 
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Fig. 1. Stokes profiles along the symmetry line connecting spot 
center and disc center. The limb side is down, the center side up. 
Bottom row. intensity, /. Middle row. Total linear polarization, L. 
Top row. Circular polarization, V. Left to right: IR spectral lines 
at 1564.8 nm, 1565.2 nm, VIS lines at 630.15 nm, 630.25 nm, 
and Ti I at 630.37 nm. 



flow, whereas in the space in-between the flow velocity is 
strongly reduced dTritschler et alj 120041 ; lLanghans et alj 120051 ; 
Rimmele & Marino 2006). Different, and in some cases, con- 
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Fig. 2. Examples of Stokes V profiles in the neutral line. Top row: Stokes V of the Fe I line at 630. 15 nm. Bottom row. The co-spatial 
profiles of Fe I at 1564.8 nm. The red line gives the corresponding best-fit profile of the inversion. The last column shows a profile 
in the center side penumbra. 



tradictory results were found on the correlation between these 
flow filaments and intensity filaments (cf. the discussion in 
lLanghans et"aT]|2005l) . 

Using spectropolarimetric observations of sunspots, the in- 
formation content of the observations was extended to encom- 
pass both the magnetic field topology and the flow fields at the 
same time. The problem of the penumbral fine structure how- 
ever was not resolved by this, but got more complex again. 
The Evershed flow could be shown to be close to horizon- 
tal, but the average ma gnetic field turned out not to be (e.g. 
Bello t Rubio et al] [2003). A possible explanation for the mis- 
match between flow and field direction in an ionized plasma 
is the assumption that at the spatial resolution of around 1", 
achieved in current spectropolarimetric observations, two differ- 
ent magnetic field components are seen inside the same pixel. 
These two magnetic field components have a different topology 
and flow field each. 

One way to visualize this is to take a cut along the symme- 
try line of a sunspot, the line that passes from disc cen ter across 
the spot center towards the clo sest solar limb position dColladosI 
12003: iBellot Rubio et al.ll2004l) . Figure Q] shows profiles of in- 
tensity, total linear polarization, L = V2 2 + U 2 , and circular 
polarization (Stokes V) along the symmetry line for one of the 
observations analyzed in the present study. Both circular and lin- 
ear polarization show a neutral line on the limb side and center 
side, respectively. At the location of these neutral lines, the pro- 
file shape changes abruptly. On the limb side across the neutral 
line of Stokes V, stronger fields with smaller line shifts turn into 
weaker fields with larger line shifts. On the center side across the 
neutral line of L, the same change from strong fields with slower 
flows to weaker fields with stronger flows is seen. 

A more direct evidence of these multiple magnetic field com- 
ponents is given by the Stokes V profiles close to or in the neu- 
tral line of V (Fig. |2]l on the limb side. The neutral line indi- 
cates a local minimum of the Stokes V signal amplitude, because 
most of the field lines are perpendicular to the line of sight. In 
magnetograms, the retrieved signal passes through zero in the 
neutral line, because the polarity of the magnetic field changes 
from, e.g. positive to negative, but in spectra the V signal does 
not disappear completely due to the multi-component penumbral 



structure (cf. ISchlichenmaier & Collados 2002). These profiles 
in the neutral line show not only two cr-components in the cir- 
cular polarization signal, but several local extrema (both minima 
and maxima). To generate such complex patterns, the magnetic 
field has to be complex as well. Appendix [A] shows how such 
complex profiles can be created by the addition of two regular 
Stokes V profiles with opposite polarities. This simple addition 
of profiles however does not tell, if the two magnetic field com- 
ponents leading to the V profiles are actually located in the same 
pixel. If the spatial resolution of the polarimetric observations 
is insufficient to separate the flow filaments seen in high res- 
olution observations from their surroundings, field components 
from two different spatial locations would be added up. 

There however is an indicator that suggests that the field 
components are not spatially separated in the horizontal dimen- 
sions, but vertically interlaced, the Net Circular Polarization 
(NCP). A non-zero value of NCP - as observed in the 
penumbra of sunspots (e.g. ISanchez Almeida & Litesl 1 992; 
Sola nki & Montavonl 11993b iMartmez Pilletl 120001: iMiiller et all 
l2002h iMiiller et alJ 12006) - indicates discontinuities or at least 
strong gradients of field properties along the line of sight inside 
the formation height of spectral lines. 

Most of the observed properties of profiles in the penum- 
bra of sunspots can be explained with the uncom bed penumbra 
model suggested bv lSolanki & Montavonl d 1 993b - In this model, 
a more vertical field component winds around horizontal flow 
channels, which carry the Evershed flow. The success of this 
model is related to the fact that it is able to explain the most 
prominent peculiarities of sunspots: a horizontal flow field in a 
non-horizontal field topology, the appearances of neutral lines in 
total linear polarization, L, and Stokes V, and the NCP. However, 
the model itself gives no reason, why its spe cific topology should 
exist i n the penumbra. The simulations of ISchlichenmaier et al.l 
dl998l) showed that such a configuration of embedded flow chan- 
nels results from the temporal evolution of flux initially located 
at the bounda ry layer between the su nspot and the field-free 
surroundings. Thomas & Weissl ((2004) suggested that turbulent 
convection outside the spot pulls down field lines and thus pro- 
duces the filamentary structure of the sunspot. In both cases, 
the penumbra consists of simila r embedded flow channel s, but 
created by different mechanisms. Spru it & Scharmerl d2006l) and 
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Table 1. Spectral lines observed with POLIS (top) and TIP (bot- 
tom). 



Fig. 3. Continuum intensity maps of NOAA 10425 in the near 
infrared on 7.8.2003 (left) and 9.8.2003 (right). The distance be- 
tween tick marks is 1 arcsec. 



IScharmer & SpruiH (120061) recently suggested a completely dif- 
ferent model, where the intensity filaments are related to the ex- 
istence of field-free gaps below the visible surface layer. 

In this paper, I want to show that a rather simple uncombed 
model using two depth-independent magnetic components is 
sufficient to reproduce simultaneous observations in infrared and 
visible spectral lines with their different responses to magnetic 
fields. The inversion setup is described after the observations and 
data reduction (Sect.|2]l. I study average field properties and their 
radial variation in Sect. [4] I then deduce a geometrical model 
of the sunspot by integrating the surface inclination of the fields 
(Sect.[5]l. The physical properties of the model like field strength, 
velocity, or field orientation, correspond to the best-fit model at- 
mospheres necessary to reproduce the observed spectra. I shortly 
investigate the temporal evolution in Sect. [6] with emphasis on 
which parts change with time and which do not. In the discus- 
sion (Sect. |7]l, I address the question of penumbral heat transport 
in the context of the results of the previous analysis. AppendixlAl 
shows how complex profiles can be created with simple assump- 
tions. Appendix FBI shows an example of the integration along a 
single cut through the penumbra and the full 3-D model from 
different viewing angles. Appendix ICl describes the analysis of 
a time series of speckle-reconstructed G-band images used to 
derive the characteristic time scale of penumbral intensity varia- 
tions. 

I used a simple inversion setup with field properties constant 
with optical depth, corresponding to a horizontal separation of 
the field components. A more sophisticated setup, taking into 
account the variations along the line-of-sight, is planned to be 
discussed in a forthcoming paper. 

2. Observations & Data reduction 

For this work, two observations of the sunspot NOAA 10425 
were analyzed. They were taken on the 7th and 9th of 
August 2003 with the two vector spectropolarimeters of the 
German Vacuum Tower Telescope (VTT) on Tenerife: the 
Polarimetric Littrow Spectrograph (POLIS : iBeck et al.ll20 05b) 
and th e Tenerife Infrared Polarimeter (TIP; Martinez Pille t et al.l 
1999). The two instruments were used simultaneously using an 
achromatic 50-50-beamsplitter. During the observations, POLIS 
was remote-controlled by TIP to ensure strictly simultaneous ex- 
posures. The image motion was reduced using the Correlation 
Track er System (Schm idt & Kentischerlll995H Ballesteros et all 
119961) . The Stokes vector in the visible and infrared spectral lines 
listed in Table [TJ was retrieved. An integration time of around 
3.5 sec was used for each slit position. The slit width was 0736 
for TIP and 075 for POLIS. The scanning step width was 07 36. 
Spatial sampling along the slit was 0737 for TIP and 07145 for 
POLIS. The POLIS data were interpolated later to have the same 
spatial sampling along the slit as TIP. Figure|3]displays maps of 
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(c) lBeflot Rubio et al. 1 2000); (d) L. R. Bellot Rubio, priv. comm. 
All values are given as used in the inversion; they may slightly deviate 
from the sources cited in A and log(g/) due to changes in the adopted 
solar iron abundance and new IR line measuerements. 



the continuum intensity in the infrared. The 2-D intensity maps 
were constructed from the intensity along the slit for each scan 
step, the scanning direction is left to right. The heliocentric an- 
gle of the spot on the August 7th and 9th was 7° and 30°, re- 
spectively. The spatial resolution was estimated from the spa- 
tial Fourier power spectrum to be around 1". The data from TIP 
and POLIS were treated with the flatfielding procedures and po- 
larimetric corrections for instrumental effects (see for example 
Bec k et al. l l2005bllaT) . Residual crosstalk between the different 
Stokes parameters was estimated to be on the order of 10~ 3 / c . 
The remaining rms noise in the profiles in continuum windows 
was 4 x 10~ 4 / c for the IR spectra and 10~ 3 / c for the visible spec- 
tra. 

The dat a alignment was done in the sam e way as describe d 
in detail in lBeckl d2006l) or the Appendix of IBeck et all d2007l) . 
The wavelength scale was also set up like in the latter, with 
the blue-shift values predicted by the quiet Sun (QS) model of 
Borr ero & Bellot Rubiol (12002b as a reference. The Stokes pro- 
files of each spectral range were normalized to the continuum 
intensity of the quiet Sun at disk center in a two-step procedure. 
An average QS profile for each wavelength range was calculated 
from pixels located outside the sunspot. This profile was nor- 
malized to unity with its respective average continuum intensity. 
The off-center position was taken into account by a multiplica- 
tion of the QS profile and all other profiles with the appropriate 
limb-darkening coefficient. 



3. Inversion of the Stokes profiles 

The four visible and the three infrared lines have been in- 
verted together using the SIR code (Stokes Inversion based 
on Response fu nctions; iRuiz Co bo & del Toro Iniestal 119921; 
Ruiz Cobo1 ll998l) . To facilitate a good choice for the inversion 
setup and the initial model components to be used, I created 
masks of the inversion type. I used the intensity maps in infrared 
to define the boundaries between umbra and penumbra, and 
penumbra and surrounding granulation, respectively. Outside the 
sunspot, I used a threshold in polarization degree to distinguish 
between field-free pixels and those with a polarization signal 
sufficiently large to derive the magnetic field. The threshold was 
set to 0.4 % for the infrared spectral lines (1564.8, 1565.2 nm) 
and 0.75 % for the visible spectral lines (630.15, 630.25 nm). If 
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Fig. 4. Comparison of the observed (+) and best-fit profiles (solid line) for one spatial position in the neutral line of Stokes V. Left 
to right: infrared lines at 1564.8 nm and at 1565.2 nm; visible lines at 630.15 nm, 630.25 nm and 630.37 nm. Top to bottom: Stokes 
IQUV as fraction of the continuum intensity. The black vertical lines mark the zero-wavelength of the respective wavelength scales. 
The dispersion on the x-axis is in mA for all spectra. 



any one of the lines exceeded its respective threshold, a magnetic 
field was assumed to be present. 

All model atmospheres were prescribed as a function of con- 
tinuum optical depth, r, in the range from log r = 1 to -4. For 
the initi al temperature stratific ation, I always adopted the HSRA 
model (Ginge rich et al.|[l97ll) . Temperature was allowed to be 
varied with two nodes, i.e. perturbations of the initial model at- 
mosphere with a straight line of arbitrary slope were possible. A 
contribution of stray light to the observed profiles was also al- 
ways allowed for; as proxy of the stray light profile the average 
QS profile was used. 

I employed a two-component inversion with one magnetic 
and one field-free component for pixels outside the sunspot with 
significant polarization signal. The free parameters of the field- 
free component were the temperature, T, and line-of-sight (LOS) 
velocity, v. For the magnetic atmosphere component, additional 
parameters were the field strength, B, the LOS inclination, y, and 
the azimuth of the field in the plane perpendicular to the LOS, 
if/. Except for temperature, all atmospheric parameters were as- 
sumed to be constant with optical depth. Those spectra, whose 
polarization degree was below the threshold, were not inverted. I 
assumed a single magnetic component in the umbra. For all pix- 
els in the penumbra, the inversion setup used two independent 
magnetic components in each pixel. The macroturbulent veloc- 
ity was the only parameter that was forced to be equal in both 
components. All quantities besides T were again constant with 
depth. 

For SIR, the two inversion components are fully equivalent; 
the changes applied to their initial model atmospheres are only 
driven by the need to minimize the deviation between observed 
and synthetic profiles. In the inversion of the spectra of neighbor- 
ing pixels the roles of the two inversion components thus may 
be exchanged, i.e. component 1 may show a large flow veloc- 
ity and component 2 is at rest in one case, whereas on the next 
pixel it is the opposite way. The results of the inversion have 
thus to be sorted somehow to provide a smooth spatial variation. 
I used the inclination to the surface as criterion to separate the 
two inversion components into different maps: the more vertical 
component of each pixel is assumed to be the static background 
(bg) field of the spot, and the more inclined component to be the 
flow channels (fc). This selection by a single criterion can lead to 
ambiguities, when the inclination of the two components is simi- 
lar. Near the umbra-penumbra boundary, the field inclinations of 
the two components were nearly identical for both observations 



of the sunspot (cf. Beck 2006, Fig. 5.14). In this area, I modi- 
fied the criterion and used the temperature of the components: 
the hotter component was chosen to be the bg component, the 
cooler one the fc. This yielded smooth temperature maps, which 
otherwise showed clear indications that the identification of the 
components by inclination was wrong. 

Figure |4] shows an example of observed and best-fit profiles 
for one location in the neutral line of Stokes V. The weak Fe I 
line at 630.35 nm has been left out in the graphs to save some 
space. The inversion is not able to reproduce asymmetric pro- 
files, and thus fails to retrieve the observed profiles down to 
the finest details (e.g., Stokes Q and U of the visible lines). It 
however still catches the overall shape of the profiles quite well. 
Especially in Stokes V one and the same atmosphere leads to 
strongly differing profiles in, e.g., 1564.8 nm and 630.15 nm. 
The good agreement of observed and best-fit profiles in both 
wavelength ranges using constant field properties is related to the 
only small differ ence in formation height be tween the infrared 
and visible lines. Cabre raS plana et alj d2005l) studied the prop- 
erties of several photospheric lines and concluded that the forma- 
tion height difference between the 1 .56 fi lines and those at 630 
nm is smaller than 100 km. The profiles shown in Fig.|4]have a 
non-zero NCP; the NCP is not reproduced by the best-fit profiles. 
Even if the observed sunspot has a non-zero NC P with differ- 
ent b ehavior in infrared and visible spectral lines dMiiller et al.l 
2006), the main information contained in the spectra is at first the 
average properties of the magnetic field like field strength. The 
NCP then yields information on the v ariation alon g the LOS, but 
always around the average value. In iBeckl d2006l) I found that 
the analysis of the same data set using an uncombed inversion 
model yielded almost identical average field properties, i.e. field 
strength or field orientation were not influenced by the choice 
of the inversion setup with or without reproduction of the NCP. 
This does not come surprisingly, because the inversion with con- 
stant magnetic field properties already reproduces the observed 
(complex) spectra fairly well (cf. Figs. [2] and 0). 

4. Results 

Prior to further analysis, the inversion results were transformed 
from the LOS reference frame into the local reference frame 
(LRF). The LRF is defined such that z corresponds to the surface 
normal and increases with height, while the x-axis points from 
the center of the spot towards the disc center. The 180° azimuth 
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Fig. 5. Top row, left to right: magnetic field azimuth of the background component (bg), field azimuth of the flow channel component 
(fc), field strength of bg, field strength of fc, and difference of the field strengths (bg - fc). Middle row, left to right: field inclination 
of bg and fc, temperature of bg and fc, difference of temperatures (bg - fc). Bottom row, left to right: LOS velocity of bg and fc, stray 
light contribution, filling fraction of bg and fc. Inclination and azimuth are given in the LRF frame; velocities are in kms~' with 
negative velocities pointing towards the observer. The velocity limits in parentheses refer to the flow channel component. White 
contour lines denote the inner and outer boundary of the penumbra. The black arrow in the azimuth map points towards disc center. 



ambiguity was resolved by assuming a radial field orientation in- 
side the sunspot, and thus always choosing the azimuth solution 
closer to radial orientation. 

4. 1 . 2-D maps of field parameters 

Figure [5] shows the resulting field azimuth and inclination to 
the surface normal. Only the results of the observation on the 
9th of August at 30° are shown here and in the following, as 
the other observation yielded similar results. The two inversion 
components have been separated by the criteria described in the 
previous section. There is an easy way to control if the selec- 
tion is reasonable. The component selected as "flow" channels 
should also exhibit significant flow velocities. And indeed this 
is the case: the flow channel component shows LOS velocities 
up to some kms~', whereas the background component is al- 
most at rest (cf. bottom row of Fig . [5j ■ This gives confidence that 
the other quantities from the inversion can really be ascribed to 
"flow channels" and "background component". 

Figure [5] also shows that the background component is 
stronger almost throughout the whole penumbra, whereas at 
the outer boundary the field strength becomes comparable. The 
patch of increased field strength in the fc component (red color) 
below and right of the umbra near its boundary to the penumbra 
is co-spatial to the locations, where the field inclinations of the 
two components were nearly identical. The strong difference to 
the field strength further radially outwards could indicate that in 
the spectra of these locations no clear signatures of two magnetic 
components were present. The temperature maps (at log r = 0) 
in the middle row show the background component to be hotter 
than the flow channels all throughout the penumbra. The stray 
light contribution to the best-fit profiles increases from 10 % 
in the umbra to 15 % in the penumbra, and strongly increases 
at the outer spot boundary. Inside the penumbra, flow channels 
and background component have a comparable filling fraction of 
around 50 %. 

Radially oriented structures ("spines", iLites et aill 1993b can 
be seen on the center side in the background components' field 



strength, the flow channel temperature, and in the relative fill- 
ing fraction of the background component. The spines start with 
locally enhanced field strength of the background component 
in the umbra, and maintain stronger and less inclined fields 
throughout the penumbra. The filling fraction of the flow chan- 
nel component inside the spines is only 50 %, compared to about 
70 % outside the spines. This indicates that the Evershed motion 
may be (at least partially) suppressed in the spine areas. Their 
continuum intensity is reduced mainly in the inner penumbra; in 
the outer part the spines are less prominent in intensity. 

4.2. Azimuthal averages 

In order to investigate the mean properties of the two compo- 
nents, I averaged the field parameters azimuthally over ellipses 
with increasing radius. In the umbra the results of the one- 
component fit are used for the background component. In Fig. [6] 
the results for the azimuthal averages of both observations (7th 
and 9th) are shown for comparison. Starting in the lower right, 
the intensity in visible and infrared continuum wavelengths has 
been used to place the transition from umbra to penumbra at 
r/r S pot—QA, where r spot is the spot radius of around 13 Mm. The 
contrast of umbra and penumbra is larger for the visible wave- 
length range. The middle right panel shows that the relative fill- 
ing fraction of bg field and flow channels changes in the mid 
penumbra; the contribution of flow channels increases from 50 
% at r/> spot =0.65 to a maximum of around 70 % at r/r spot =0.8. 
In the azimuthally averaged absolute LOS velocity the observa- 
tion at 30° heliocentric angle shows a large difference of 2 kms 
between bg component and flow channel, whereas at 6° the sig- 
nature of the Evershed flow is negligible. In both cases the bg 
component shows increasing LOS velocities towards the outer 
boundary. 

The field strength and field inclination of both observations 
are similar, indicating little change of the sunspot field topology 
in two days. In both observations, the bg component is stronger 
by 0.5 kG in the inner and middle penumbra, whereas at the outer 
boundary the strength of both components is nearly identical. 
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Fig. 6. Le/f column, top to bottom: Radial variation of B, LRF in- 
clination, temperature at log r = 0. SoZ/of: bg component, dashed: 
flow channel component. Right column, top to bottom: averaged 
absolute LOS velocity; stray light amount (solid) with overplot- 
ted filling fraction of the fc component (dashed); intensity of vis- 
ible (solid) and infrared continuum (dashed). In this case, black 
(grey) lines refer to the observation on 9th (7th) of August. The 
x-axis gives the fractional radial distance r/r spot , the black ver- 
tical lines mark the inner and outer boundary of the penumbra. 
The horizontal black line in the inclination plot marks 90°. 



The slope with radius of the flow channel field strength decreases 
in the mid penumbra at r/r spot = 0.7 (dotted lines), but the field 
strength continues to drop with radius in both bg component and 
flow channels. The inclination of the flow channel component on 
average exceeds 90° for r/r spot > 0.9, whereas the bg component 
never turns into horizontal fields. Its maximum inclination at the 
outer penumbral boundary is close to 60°. 

The temperature plot in the lower left panel shows that the 
radial variation of the temperature in the background component 
follows closely the radial curve of intensity. The temperature 
curve of the flow channel component has a similar shape with 
reduced amplitude, but is displaced towards the outer penum- 
bral boundary relative to the intensity or bg temperature curve. 
The decrease of the temperature in both components at the outer 
penumbral boundary is presumably connected to a trade-off be- 
tween stray light and temperature in the inversion. At the outer 
boundary, the intensity level and shape of the profiles allows to 
use a larger stray light amount to reproduce the observed spec- 
tra, whereas in the umbra and penumbra the QS profile simply 
does not fit to the s pectra. These results are in good agreement 
with th e findings of iBorrero et al.l (l2004t) or lBellot Rubio et al.1 
(120041) . both in the absolute values of, e.g., field strength or field 
inclination, and in their radial variation. 



4.3. Field-aligned flows 

With the assumption that on large spatial scales the flow 
velocities in the penumbra are only due to the radially 
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Fig. 7. Top: Horizontal and vertical velocities for both inver- 
sion components. Black: bg component. Grey: flow channel 
component. Dashed: vertical velocity. Solid: horizontal velocity. 
Bottom: field inclination, y (solid), and flow angle, a (dashed). 



aligned Evershed flow, the LOS velocity can be decom- 
pose d into its horizontal, Vh(r), and vertical component, v 7 (r) 
(e.g.lSchlichenmaier & Schrmdtll20"oot iBellot Rubio et al.ll2003l 
2004). This allows to derive the flow angle, a(r), i.e. the incli- 
nation of the flow direction to the surface, at a given radius. The 
flow angle can be derived separately for each inversion compo- 
nent, and then be directly compared to the average field inclina- 
tion of the component, y(r). Figure|7]displays the corresponding 
results for the observation on 9th of August. The flow channel 
component shows upfiows of around 1 kms -1 (v z > 0) in the in- 
nermost penumbra, which turn into nearly horizontal flows with 
constant 4.6 kms~' all throughout the penumbra, and finally ex- 
hibit a small downflow component for r/r spot > 0.9. The back- 
ground component shows a significant LOS velocity near the 
outer penumbral boundary, which could also be already se en in 
the azimuthal averages of Fig. [6] Bellot Rubio et al. (2004) have 
pointed out that this non-zero velocity of the background compo- 
nent is a necessary ingredient for correctly predicting the sign of 
the net circular polarization in the center side. The lower panel 
of Fig. [7] displays that in the outer penumbra the flows in the 
flow channels are parallel to the field inclination to a high de- 
gree. In the inner penumbra and for the background component 
the agreement is worse, but this may also be due to the intrin- 
sic inaccuracies of the method. The azimuthal fine structure of 
dark and bright filaments is ignored by assuming that the abso- 
lute flow velocity v = v(r) only depends on radius. Given these 
limitations, I think that the result supports the concept of field- 
aligned flows: upfiows (downflows) are seen where the magnetic 
field inclination is below (above) 90°. 

4.4. Signature of hot upfiows in the mid penumbra 

To investigate the radial structure of the fields without spatial 
averaging, I took the values of several quantities along a single 
column of the 2-D maps at around 0° field azimuth. Figure [8] 
shows this column and its surroundings. At about the middle of 
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Fig. 9. Parameters along the column marked in Fig. [8] 
Left column, top to bottom: continuum intensity, tem- 
perature, field strength. Right column, top to bottom: 
field azimuth, field inclination, line-of-sight velocity. In 
all graphs besides intensity, solid lines indicate quanti- 
ties of the background field and dotted lines those of the 
flow channel component. For intensity, solid and dotted 
correspond to visible, respectively, infrared continuum 
intensity. The line-of-sight velocity of Ti I is given by 
the long-dashed line in the velocity graph. The vertical 
black line marks the outer penumbral boundary. The 
second hot upstream intersected by the cut at d=6.2 
Mm has a similar signature in all quantities as the up- 
flows at the inner penumbral boundary. 




Fig. 8. (Top row): continuum intensity and temperature of the 
flow channel component. (Bottom row): velocity and inclination 
of the flow channel component. The white lines mark the col- 
umn, along which the parameters shown in Fig. [9] were taken. 
Black arrows point towards hot upstreams in the mid penumbra. 



the penumbra a hot upstream is located. It has a strong signature 
in most quantities along the cut shown in Fig. [9] The line-of- 
sight velocity in the flow channel component drops from 2 kms -1 
further inwards to slightly below zero; it takes around 2 Mm to 
speed up again. The same is seen in the line core velocity of 
Til. The inclination of the field is reduced from 88° t o 55°, but 
j umps to a round 77° on the next pixel ( = (X'37) already. Rimmele 
(2004) and Rimmele & Marino (2006) also found a rapid change 
of upflows into horizontal outflows on scales of around Qf.'6. The 
temperature of the flow channel component increases by around 
1000 K on the upstream location. In the continuum intensity, an 
increase is seen near the upstream, but displaced outwards by 1 
Mm. Several similar hot upstreams in the mid penumbra with a 
change of temperature, flow velocity and field inclination can be 
seen in Fig. [HJ two examples near the left and right border of the 
section shown are highlighted by arrows. 

Significance of variation The inversion procedure treats the 
spectra of each position independently of their surroundings; 
the inversion components are then separated with the two cri- 



teria (inclination, temperature) as described in Sect. [3] This can 
produce some strong pixel-to-pixel variations in the 2-D maps 
of the inversion parameters. I think that for the case shown in 
Fig- El I can exclude such an origin of the jump in atmospheric 
properties. The field inclinations differ significantly; the contin- 
uum intensities and the line core LOS velocity of the Ti line 
are quantities that are not related to the inversion procedure, and 
they both show also the signature of a hot upstream: a local in- 
tensity maximum with a vanishing LOS velocity. Note that the 
situation is different fr om the "sea serpent" shape suggested by 
ISchlichenmaierl d2002l) . as I have no indications for downflows 
or field lines pointing down near the upstream. A possible con- 
figuration of the field lines using the integration of the inclination 
(cf. the next section) is shown in AppendixlBl Fig. IB. 11 



5. Derivation of the 3-D field topology 

At first sight, the inversion with two magnetic components, 
whose properties are constant with optical depth, does not con- 
tain information on the vertical structure of the magnetic fields. 
Taking into account that the information on the field vector is 
available for an extended spatial region allows to derive a 3-D 
model nonetheless by integrating th e field inclination. The same 
metho d was already proposed in ISchlichenmaier & Schmidtl 
(2000), and applied to the radial variation of the flow angle de- 
rived f rom spectroscopy. It ha s also b een used in ISolanki et al.l 
( 2003;j) and Iwlegelmann et al.l (120051) on measurements of the 
magnetic field vector. ISolanki et al. " 32003b used the requirement 
that the field azimuth did not change along the integration direc- 
tion, which is fulfilled for the sunspot observed. 



5.1. Integration of field inclination 

For each radial position, r,, the azimuthally averaged inclina- 
tion, y(r,), of background component and flux tube component 
are known (cf. Fig. [6j. With the assumption that the inclination 
does not strongly change with depth below the surface (or height 
above it), the geometrical height, h, of a field line in the solar 
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Fig. 10. Integrated inclination of the data of 7th (left) and 9th 
(right) of August. (Lower thick black line): bg component; (up- 
per thick black line): fc component. (Dark grey shaded): same 
curve including inclination variations. (Thick dashed lines): the 
boundary layers of the spot model of JS94. Thin black line: 
radial variation of the infrared continuum intensity. (Vertical 
dashed lines): inner and outer penumbral boundary. The thin 
black line at around z=0 km marks the location of the r = 1 
level. The light grey shaded area at (z=0 km,r = 3 Mm) de- 
notes the formation height(s) of the spectral lines. The dotted 
lines are identical to the bg component curve, but shifted in z to 
have the inclination observed at the r = 1 level. 



atmosphere can then be calculated from an integration of the in- 
clination by 



h(n) 

■i— ' tan "v 



urn T (r,) 



Ar, 



(1) 



where Ar = 330 km is the difference of the respective semi- 
major axes of two subsequent ellipses. 

Two additional steps were applied to the curves of h(r,) for 
background field and flow channels. The flow channels are not 
observed in the umbra; hence, no inclination values are avail- 
able there. To have a common reference height in both integrated 
curves, I forced the height of bg and fc at rj r spot = 1 to be identi- 
cal. This was motivated by the fact that spectra at the outer white- 
light boundary still showed the signatures of multiple compo- 
nents. Thus, both flow channels and background component 
have to be present within the formation height of the spectral 
lines, and hence, must be present at the same geometrical height 
as well. I also assumed that the isosurface of T500 = 1 has a slope 
of 3 deg w hen going from the inner towards t he outer penumbral 
boundary ( Schlichenma ier & SchmidJl2000h . corresponding to 
a Wilson depression of 380 km in the um bra. The Wilson de - 
pression agrees with the number given by Math ew et alj ([2004), 
who however found no smooth radial change, but a pronounced 
jump of 280 km at the umbral-penumbral boundary. Figure [TOl 
displays the finally resulting curves for both observations on 7th 
and 9th of August. Both observations yield very similar results. 

The integrated inclination curve of the flow channel com- 
ponent yields a slightly elevated arched loop. The apex height 
above the r = 1 level is around 300-400 km. With the height 
fixed to km at the spot radius, the curve intersects the r = 1 
level at around r/r spor = 0.6. Outside the sunspot, the field 
lines point downwards. As the background component incli- 
nation never reaches horizontal fields, the integration yields a 
much steeper curve than for the flow channels. If taken at face 
value, the curve reaches a depth of 2 (4) Mm at r/r spor = 
0.8(0.6), implying a thick penumbra, not a shallow surface layer. 
For the background component fortunately a direct comparison 
with a theoretical model is possible. I overplotted the bound- 
ary layers between umbra and penumbra and between penum- 
bra and surroundings fr om the magneto-static sunspot model of 
Uahn & Schmidl dl994L JS94) in Fig. [lOj I reduced the radius 
in their original calculation to 93 % and 95 %, respectively, of 
its value, to fit with the dimensions of the sunspot in the ob- 
servations, and shifted the curves in height to be at z=0 km at 
r I rspot - 1 like the integrated inclination curves. For rjr spot be- 
tween 0.6 and 1, the agreement between the integrated bg curve 
and the magneto-statical model is astonishingly good, whereas 



for smaller radii the JS94 model is steeper than the integrated 
curve. 

The good agreement in the mid to outer penumbra comes 
a bit as a surprise: the JS94 model gives the location (and thus 
also the field inclination) of the boundary layer between spot and 
surroundings at depths up to some Mm, whereas the integrated 
curve uses the inclination observed at the surface close to the 
t = 1 level. The observed surface inclination thus seems to be 
identical to the inclination in the deeper layers - which was sim- 
ply an assumption in the derivation of the integrated curves. The 
deviation between the curves also points in the correct direction: 
if the inclination changes with depth, the fields should get more 
vertical in the deeper layers, as the expansion of flux concentra- 
tions happens near the surface layer. Thus, the surface inclina- 
tion should be larger than in the deep layers, leading exactly to 
the less steep integrated curve as obtained. 



Limitations of the integration Despite the surprisingly good 
agreement with a theoretical model - which also indicates that 
the main assumption of a small depth dependence of the incli- 
nation could be valid - the integration of the surface inclination 
of course can not be fully correct. In the left panel of Fig. [10] 
I overplotted the approximate formation height(s) of the ob- 
served spectral lines throughout the penumbra. All information 
retrieved from the observed spectra thus only refers to this small 
layer of the atmosphere. There is no guarantee that fields do not 
for example bend strongly as soon as they leave the height range, 
in which the spectral lines are sensitive. In fact, one expects ex- 
actly this behavior for the field lines that pass the upper bound- 
ary of the formation hei ght, because of the exponentia l decrease 
of density with height. Sanc hez Cuberes et alj ([2005) found no 
height dependence of the field inclination using spectral lines 
similar to the 630. 15/25 nm pair, but this would be valid only in- 
side approximately the same formation height range as plotted in 
Fig. [10] To investigate the influence of inclination changes with 
height (or depth) on the integration of inclination, I repeated the 
integration with the assumptions that the inclination may vary 
with height by +20 % for the background component, respec- 
tively, + 15 % for the flow channels. The upper and lower limits 
of the retrieved curves are given by shaded areas in Fig. [10] As 
the expected change should be a deviation towards more inclined 
fields in higher and less inclined fields in lower layers, any curve 
that would not leave the shaded areas could be in agreement with 
the observed surface inclinations, within the 15-20 % range of 
variation allowed for. Even if one takes the extreme case of a 
field line starting at the lower and ending at the upper boundary 
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Fig. 11. Top left: Image of NOAA 10425 in the G-band from the 
DOT on La Palma, taken on 9th of August about half an hour 
after the observations at the VTT. The white arrow points to- 
wards disc center. The white rectangles mark regions with an 
identically structured umbra-penumbra boundary in both inten- 
sity maps. Top right: infrared continuum map. Bottom: Top view 
of the 3-D model. See text for the description. 



of the shaded area, the curves would not change enough to con- 
tradict the description above. The background component still 
would imply a thick penumbra; the flow channels could change 
from non-elevated to strongly elevated loops, but without giv- 
ing rise to a new topology of the flow channels relative to the 
background component. 

5.2. A three-dimensional model 

In the previous section, the integration was performed for az- 
imuthal averages, where the azimuthal fine-structure of the 
penumbra is lost. The integration can however also be performed 
for individual radial cuts. To simplify the calculation and to 
smooth out slightly the pixel-to-pixel variations in the inversion 
results, I decided to use 92 bin^j of about 4 deg angular extend 
for the construction of a 3-D model that takes the azimuthal fine- 
structure into account. The model itself is represented in the fol- 
lowing way: the integration of the background component gives 
a mesh of height with radial and azimuthal position, /z(r,, <pi). 
This mesh is interpolated to a smooth surface, whose color is 
set to represent the field strength of the background component. 
The flow channels are overplotted as thin lines, with a color code 
corresponding to their temperature. The temperature is scaled in- 
dividually between the respective minimum and maximum for 
each of the 92 bins; thus, the color bar only gives an average 
range of temperature. This description applies to Figs. [TT] and 

fiO 

As an intermediate step to the full 3-D model, Fig.fTTIshows 
a top view of the model. This view corresponds to the com- 
monly used 2-D maps of physical quantities. For comparison, 
I also show a speckle-reconstructed image of NOAA 10425 in 
the G-band taken with the Dutch Open Telescope (DOT) on 9th 
of August 2003 about half an hour later than the observation 
analyzed here. The figure serves as a comparison of the spatial 
resolution of the polarimetric to speckle-reconstructed data. The 

1 It needed to be a multiple of 4 for technical reasons. 
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Fi". 12. The temporal evolution of the penumbra with 7 min ca- 
dence. Data of NOAA 10425 on 9th of August 2003, 9:36 to 
10:40 UT. The lower row is the infrared continuum intensity, 
the upper row the co-temporal and co-spatial map in the G-band 
from the DOT. Tick marks are arcsec; time increases from left to 
right in each row. 

boundary between umbra and penumbra has the same shape in 
both observations (cf. inside the white rectangles), even with the 
time difference of half an hour. Bright penumbral grains can be 
seen in the DOT map near the umbral boundary, but only on the 
limb side. 

Figure I5~2l shows the full 3-D model for both observations 
for three different viewing directions^. I emphasize that the 
model is derived more or less straightforward from the observed 
spectra, with in principle the single assumption that the integra- 
tion of the field inclination is reasonable. If the model properties 
are projected to the surface at z=0 km, one has exactly the vector 
field and thermodynamic properties that gave the best-fit to the 
observed Stokes profiles in - including the weak line blends - 
3 infrared and 4 visible spectral lines. The model(s) show some 
interesting peculiarities and give rise to several interesting ques- 
tions that could be addressed by analyzing them. The observa- 
tion close to disc center (7.8.2003, 6° heliocentric angle) is very 
uniform all around the spot, whereas in the other case there is 
a much larger difference between center and limb side, where 
the flow channels appear to be much more elevated. The back- 
ground component exhibits a ragged subsurface structure, where 
less inclined and steeper fields coincide with an increased field 
strength. This could be c losely related to the question if these 
spines ' dLites et alj [l993) are co-spatial to darker or brighter 
filaments. In general, a study on the relation of intensity to this 
specific geometry would be interesting, because positive as well 
negative correlations between field strength, inclination, and in- 
tensity have been found up to now. Before an intense study on 
the question of the validity of the integration of inclination has 
been performed, I refrain from speculating further at this point. 
The 3-D models sho w some sim i larity to the artists sketch of a 
sunspot in Fig. 4 of I Weiss et alj (2004), but I remark that they 
are derived (directly) from observations here. 

6. Temporal evolution 

To investigate the temporal evolution, I used a time-series taken 
after the map on 9th of August. Figure [12] shows a small sub- 
section of the full field-of-view (FOV) at the DOT together with 
a co-temporal observation taken at the VTT from 9:36 to 10:40 

2 Two animations of the model for the 9.8.2003 are available in the 
online material. 
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UT, after the map analyzed to derive the field topology. This is 
the upp er part of the FOV containing the sunspot which was not 
used in lBeck et al.l d2007l) . where the data properties and align- 
ment method are discussed in more detail. The images show 
a peculiarity of the penumbral dynamics that I think has not 
been given enough attention up to now: the global structure of 
the sunspot, i.e., the boundary between umbra and penumbra 
marked with a white rectangle , does not evolve at all in around 
1 hour, whereas the intensity pattern of bright penumbral grains 
(PGs) evolves much faster on scales below 30 minutes (cf. also 
Appendix 0. This suggests that the dynamical evolution does 
not change the geometry. In the context of the previous sections, 
this could be interpreted as a static background component with 
a dynamic flow channel component. 

7. Discussion 

The penumbra has been subject of many studies, which in most 
cases did not agree in their results. The observational findings 
are converging in some points, which I think I can substan- 
tiate further by the results of this investigation. The observa- 
tions agree that the topology of the penumbra is complex. The 
Evershed flow happens along the nearly horizontal flow chan- 
nels, whereas the averag e field inclination is not horizontal (e.g. 
iBellot Rubio et al.ll2004 . The Stokes V profiles of Figs. and g] 
clearly show that at a resolution of 1" at least two independent 
magnetic components are present in each pixel. This could be 
an artifact due to the spatial resolution, but if the picture of the 
uncombed penumbra is valid, it will be the case even if the flow 
channels were fully resolved. As long as the flow channels are 
not optically thick and located inside the formation height of a 
spectral line, both horizontal and more vertical fields would be 
seen at the same time. Stokes spectra with higher spatial resolu- 
tion from, e.g., the recently launched HINODE satellite could be 
used to see if the signature of multiple components disappears at 
high spatial resolution. 

Some results turn out to be identical regardless of the in- 
version method employed, being it a simple or complex model. 
If polarimetric profiles are analyzed in terms o f two compo- 
nents (or also one component, but with gradients dBorrero et alJ 
120041) '). one retrieves a stronger, less inclined component with 
small flow velocities, and a strongly inclined one, which har- 
bors large flows. The flow is located inside the magnetic fields, 
as the flow velocity is derived from the Doppler shifts of the 
polarization signal. The same classification into less inclined 
static field and more inclined flow channels was obtained by 
lLanghans et al .] (120051) from magnetograms. All recent inver- 
sions performed agree that the more in clined component bends 
downwards in the outer penumbra (Westendorp Plaz a et a"D 
1200 It IBellot Rubio eTall l2004t iBorrero et alJ 120041) . Finally^ 
hot upflows were found at the in ner penumbral boundary in 
also some other recent studies ( Schlichenmaier et al. 2004; 
Tritschler etafl l2004t IBorrero et alJ l2005t iRimmele & Marinol 
20061 IBellot Rubio et al.ll2006l) . 



With respect to the topology of the fields, I introduced the 
concept of the integration of the surface inclination. This ap- 
proach relies on the assumption that the field inclination does 
not, or only slightly, depend on depth or height in the solar at- 
mosphere. Even if the validity of the assumption has not been 
addressed at all in this work, the method yields reasonable re- 
sults. For the background component, I find a surprisingly good 
agreement of the integrated curve with the outer boundary layer 
between the sunspot and its surroundings in the magneto-static 
model of JS94. This is especially unexpected, as their bound- 
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Fig. 13. The effects of the inclined line of sight on observations 
of the limb or center side. Thick black: integrated LRF incli- 
nation of the 2-component inversion as in Fig. [10] flow channel 
component (upper line) and background component (lower line). 
The LOS is overplotted in purple for the center side and in or- 
ange for the limb side. 



ary line gives the field inclination in deep layers of some Mm, 
whereas I used the inclination on the surface. The integration 
of the flow channel inclination yields slightly elevated arched 
loops. 

Another peculiarity refers to the behavior of the inversion 
results at around r/r spol = 0.65. At this location, the fill factor 
of the flow channels strongly increases. Westendorp Plaz a et all 
(200l j) found a local ma ximum of intensity at r/r spot = 0.6-0.7; 
Bellot Ru bio et all d2004l) found a change of inversion results in 
the mid penumbra, a jump in the quantities of their more vertical 
component and an increasing filling fraction of the flow chan- 
nel component. Interestingly, an intersection point of the surface 
with the integrated curve is located near that radius. The inver- 
sion is performed pixel-by-pixel on only the spectra from a sin- 
gle location, whereas the integration uses the inclination results 
along a radial cut to create the curve of the field lines. As these 
two procedures are fully independent of each other, I think that 
this co-spatiality is not only pure coincidence. Taking the inte- 
grated curve at face value, the obvious explanation is that flow 



channels on average cross the surface around rjr, 
then have a stronger signature in the observed spectra. 



spot 



0.65 and 



7. 1 . Line-of-sight effects 

The 3-D models derived by the integration for different helio- 
centric angles show some deviations in their structure, even if 
azimuthally averaged values are nearly identical. The model of 
the observation near disc center is very uniform, while the other 
observation shows strong differences between limb and center 
side. As both observations were analyzed in exactly the same 
way, the difference should be contained already in the spectra. 
A second point related to LOS effects is the comparison with a 
high-resolution image of the same sunspot from the DOT tele- 
scope (cf. Fig. fTTb . Bright penumbral grains can only be seen 
on the limb side. This could be due to the simple geometrical 
effect depicted in Fig. [13] If the LOS is inclined to the surface, 
on the center side a limiting angle exists, where a structure will 
be hidden by its own continuation. The two LOS lines on the 
center side in Fig. [13] intersect the bg and ft curves at locations 
where their inclination is the same as the heliocentric angle of 
the observation. The heliocentric angle of the observations on 
9th of August was around 30°; the minimum field inclination of 
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Fig. 14. Comparison of a snapshot of the Moving Tube Model (top panel, courtesy of R.Schlichenmaier) with the flow channels' 
geometry from the observation (lower panel). The plots have been arranged to have the same scaling and position in the horizontal 
axis which gives the distance from the spot center in Mm. Both panels display the flow channel component and the background 
component with their respective field strength as color coding. The color bar gives the color coding for field strength in both graphs. 
Velocities are overplotted as arrows, whose length gives the absolute velocity. The direction of the arrows for the lower plot is the 
field inclination. 



the flow channel component on the limb side was around 25° (az- 
imuthal average around 40°, cf. Fig.|6]l. Thus, the absence of the 
penumbral grains on the center side could be easily explained by 
the assumption that they are blocked by their own continuation 
as optically thick structures. 

7.2. The Moving Tube Model 

From the observational point of view, a penumbral model of 
interlocked horizontal and more vertical fields as given by the 
picture of the uncombed penumbra is in good agreement with 
observed spectra, with their general shape down to more sub- 
tle details like the n et circular polarization of Stokes V profiles 
(Miill er et al.ll2006h . From the var ious theoretical approaches , 
the Moving Tube Model (MTM) of lSchlichenmaier et al.l (1 19981) 
is the only one that results by itself in a similar geometry, in- 
cluding a flow along the horizontal channels. The MTM used 
the magneto-static sunspot model of JS94 as input for the back- 
ground component. I find that either the integrated inclination 
curve, or, to be on the safe side, the observed inclination values 
of the background field on the surface are in good agreement 
with the boundary layer between sunspot and surroundings in 
the JS94 model. However, not only the geometry of the back- 
ground is in agreement with the observations. Figure [T4l shows 
a direct comparison of the geometry, field strength, and absolute 
velocity of a snapshot of the MTM during the final steady-state 
of the simulations, and the same result as derived from the ob- 
servations. For the location of the flow channel, I took the value 
of the integrated inclination curve; for the outer sunspot bound- 
ary, I used the integrated background inclination curve. The flow 
velocity, and the field strength in flow channel and background 
component were taken from the azimuthal averages of the inver- 
sion results. The only quantity that was defined ad hoc was the 
width of the flow channels, on which the present inversion does 
not yield any information; I used a diameter of 200 km for it 



(cf. lBeckll2006l) . I find quite some similarity between the MTM 
and the results derived from the observations. Both field topol- 
ogy and the local properties (field strength, flow velocities in the 
inner and mid penumbra) do match reasonably well. Note that 
from the theoretical side the differences between the MTM and 
the integrated curve in, e.g., the height of the flow channel (~ 
100-200 km) will lead people to argue that the curve suggested 
by the integration is inconsistent with any physics. I remind how- 
ever that it reflects the azimuthal average of, at maximum, half- 
resolved flow channels, converted to a geometric height without 
much sophistication, whereas the MTM is supposed to des cribe 
an individual resolved str ucture. |Spruit & Scharmer (2006f) and 
IScharmer & Spruitl d2006l) raised the issue of the stability of flow 
channels of assumed circular diameter inside the penumbra, and 
claimed that no stable configurations are possible. The present 
investigation using model components without depth variation 
of the field properties actually does not yield information on this 
topic, but the inversion results and the integration of the inclina- 
tion with the assumption of no depth variation would still apply 
for the case of flow sheets with a small lateral and large verti- 
cal extension. Finally, the agreement of the NCP predicted by 
the MTM and the observations is already partly discussed in 
Miiller etal.1 (120061) . A more detailed comparison of observed 
NCP, the NCP predicted from the MTM, and the NCP resulting 
from a fit of an uncombed penumbral model to the same obser- 
vations used here is planned to appear in another paper. 

7.3. The "gappy" penumbra model 

ISpruit & Scharmer! (120061) and IScharmer & Spruitl d2006l) sug- 
gested a model of field-free gaps as a possible mechanism of 
the penumbral energy transport. Even if a direct comparison of 
their model with observations is difficult, as yet no (polarization) 
spectra were presented for their model, som e of their arguments 
can be compared to the present findings. [Spruit & Scharmen 
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(2006) suggested that the amount of stray light found in in- 
versions could be related to the presence of a field-free com- 
ponent. The stray light amount in the inversion of the present 
observations increases smoothly from around 10% in the um- 
bra through the penumbra (~ 10-15 %) to 80 % at the outer 
penumbral boundary. The 2-D map of the stray light contribu- 
tion (cf. Fig. [5]l shows little spatial structure besides the radial 
increase. The stray light inside telescope and the instruments 
TIP and POLIS was estimated to be around 15 % (Rezae fet al.l 
120071: ICabrera Solan aetal 1 120071) . There thus is little indication 
for a field-free component inside of the sunspot in the inversion 
results. 

IScharmer~& Spruit (2006) argue that the inversion of polari- 
metric data has a high degree of ambiguity, as similar profiles 
can result from different atmosphere stratifications. Using the 
IR lines at 1.5 ^m, I th ink that most of the ambiguities are re- 
moved. In lBeckl (120061) and lBeck et al.l (120071) I investigated the 
influence of the spectral lines on the parameters retrieved by the 
inversion. I found that e.g. field strength is restricted by the split- 
ting of the 1564.8 nm line within a limit of around +100 G. Most 
average quantities of the field topology (field strength, field ori- 
entation) are more or less uniquely restricted by the spectra; the 
main source of error are actually the spectra themselves: spa- 
tial resolution, signal-to-noise ratio, polarimetric sensitivity, and 
the polarimetric calibration. The question that however remains 
open is the 3-D organization of the magnetic fields. The present 
inversion yields field strength and orientation inside the forma- 
tion height of the spectral lines, but is not able to differentiate 
between a vertical or horizontal interweavement of field lines. 
The approach of the integration of the field inclination assumes 
coherent structures from one pixel to the next in the radial direc- 
tion, which is in my opinion highly probable, but need not be the 
case. 

The most prominent spectral feature inside the penumbra, 
the Evershed effect, is not present in the gappy model. I remark 
that all velocities derived here (besides the line core velocity of 
Til) always refer to velocities inside magnetic fields. To create 
the multi-lobed profiles in the neutral line of Stokes V, two com- 
ponents of magnetic fields with different orientation and bulk ve- 
locities are needed. If the gappy model solves the penumbral heat 
transport problems, still an explanation for the Evershed flow is 
needed. The inclination of the bg field shows an azimuthal vari- 
ation that leads to a "gappy" structure (cf. Figs. QT| and IB. 21 ). 
However, the spatial scale of t he va riation is larger than that pre- 
dicted bv lScharmer & Spruitl d2006l) . If the integrated curves are 
taken at face value, this also happens in layers far below the sur- 
face layer of t = 1 . 

Another argument in favor of the gappy model is that the spa- 
tial resolution of the observations I used may not be sufficient to 
detect the signatures of the structuring suggested by Scharmer 
& Spruit. In the other direction, it should be possible to con- 
struct a sunspot model based on their suggestions, calculate the 
resulting spectra in the 1.5 /urn and 630 nm lines, reduce the 
spatial resolution to around 1", and then invert the spectra with 
two depth-independent magnetic components. The results for the 
bg component could be compared with the p resent inversion re- 
sults. C ontrary to the regrettable sentence of Scharmer & Spruit 
( 2006) that "the agreement with observations obtained with such 
[uncombed] models is of unquantifiable significance" , I think 
it necessary to show that their model successfully reproduces 
spectroscopic or spectropolarimetric observations to support its 
validity. 



7.4. A model for the penumbral energy transport 

The long lifetimes ( on the order of 1 hr) for filaments (e.g. 
lLanghans et al. 2005), and the lack of sub merging flow chan- 
nels led Schlichenmaier & Solanki (2003) to the conclusion 
that interchange convection by rising hot flow channels is "not 
a viable heating mechanism" for the penumbra. This c laim 
has been renewed recent ly by [Spruit & Scharmer (200tJ) and 
IScharmer & Spruitl (120061) . who criticize the "paradigm" of em- 
bedded flow channels and suggest that the penumbral fine- 
structure can be explained by a model of field-free gaps reach- 
ing almost up to the solar surface. The energy transport in their 
model is then effected by convection in the field-free plasma be- 
low the sunspot. 

On the one hand, the findings of Sect.|6]support the static be- 
havior of the sunspot fields: the shape of the umbral-penumbral 
boundary and some especially dark patches inside the penumbra 
stay the same during around 1 hour. On the other hand, the inten- 
sity pattern of, e.g., the penumbral grains is completely changed 
after less than half an hour (cf. AppendixICTl. The time scales in 
the MTM model are of a comparable order. The snapshot shown 
in Fig. [14] was taken at 1 1/2 hour after the start of the simula- 
tion, but reflects the final steady-state solution. The flow channel 
in the MTM evolves rapidly in the beginning, and spans aroun d 
half of the penumbra after 30 min (Schlichenmaier et al. 1998). 

Langhans et al. (2005) used LOS magnetograms of the cen- 
ter side penumbra. These LOS magnetograms are not suitable to 
trace the dynamic evolution of the penumbra. On the center side, 
the field of the background component is parallel to the line-of- 
sight, whereas the dynamic flow channels are strongly inclined 
to it. Thus, the life times measured there reflect the slow evo- 
lution of the background component. The same applies to the 
findings of Sect. [6] the stronger field component will dominate 
the topology, and thus, the small amount of change in the geom- 
etry seen in Fig.[T2limplies that the background component does 
evolve only slowly. 

The conclusion on the impossibility of interchange convec- 
tion due to the lack of submerging flow channels needs some 
more explanations. To be convinced by the argumentation, one 
would have to agree on the fact that the penumbra is deep (some 
Mm), that flow channels originate from flux located initially on 
the boundary layer between the sunspot and its surroundings, 
and that this flux can become buoyant by heat input from the 
fully convective surroundings outside the spo{3- The first point 
is suggested by the observations, while the latter two are mainly 
based on the MTM. If one can agree on the ingredients above, I 
believe the penumbral heat transport can be effected by hot ris- 
ing flux tubes in full agreement with the observations. If flux 
becomes buoyant at the outer sunspot boundary due to heating, 
this necessarily is a repetitive process. As soon as the hot flux 
bundle has risen from the boundary layer, new different flux will 
form the boundary layer. This new flux would come from - in the 
terminology used throughout this paper - the background com- 
ponent. After a time span on the order of 30 min it would also 
have to become buoyant, and follow the previously risen flux up- 
wards. Due to the depth of the penumbra, several flow channels 
could be stacked on top of each other at the same time. The in- 
version results along a single column (cf. Sect. |4.4| and Appendix 
IB. 11 ) suggest exactly this configuration: two flow channels are 
seen along a radial cut at the same time at different locations in 
the penumbra. 



3 I guess, less than half of the solar physicists will agree on that. 



C. Beck: A 3-D sunspot model 



13 




Blocked 
convection 



Fig. 15. Schematical model of the sunspot structure. 

The question of penumbral heating then changes to the 
question if the penumbral energy losses can be replenished by 
more than one hot flow channel with a ch aracteristic repetition 
time a roun d 30 min. It has been shown by Schlichenmaier ^taTI 
(1999) and lSchlichenmaier & Solankl (f2003) that the final state 
of the MTM is not able to supply the penumbral heat require- 
ments. The final state of the MTM is however also not in agree- 
ment with the temporal evolution of penumbral fine-structure. 
Penumbral grains do not persist for hours, but fade away after 
reaching the umbral boundary, or turn into umbr al dots withou t 
a visible connection to the penumbra (e.g. Sobo tka et alJll995l) . 
The MTM never reaches this point, probably due to its boundary 
condition, which fix the lower foot point of the flow channel to 
the outer boundary layer between spot and surroundings. If the 
source of heat input is the convective surroundings outside the 
spot, it is easily conceivable that a flow channel would start to 
disappear, if it detaches completely from the boundary layer. 

The repeated ascent of flow channels of course would on the 
long term move all flux away from the boundary layer, lead- 
ing to the disappearance of the penumbra after a short time. To 
replace the flux, no submerging flow channels are needed. The 
replacement of the boundary layer can be effected by a gradual 
re-arrangement of the flux of the background component. The 
open question then is if a flow channel can become part of the 
background component, after reaching the umbra. This should 
not be impossible: the inclination difference between flow chan- 
nels and background component in the innermost penumbra is 
around 20° only. Furthermore, the inner, more vertical foot point 
of the flow channel will contribute to the magnetic field pressure 
term in the inner penumbra, pushing other field lines of the bg 
field radially outwards towards the magnetopause. 



8. Summary & Conclusions 

I have analyzed two spectropolarimetric observations of NOAA 
10425, taken two days apart, at an heliocentric angle of 6° and 
30°, respectively. The observations consisted of Stokes vector 
polarimetry of four visible spectral lines around 630 nm and 
three infrared spectral lines at 1.5/im, taken with POLIS and 
TIP. I inverted the spectra with the SIR code, using two inde- 
pendent magnetic model components. The field properties were 
assumed to be constant with optical depth. I sorted the two mag- 
netic components by their inclination to the surface normal, and 
refer to the more (less) inclined component as flow channels 
(background component). In the innermost penumbra the incli- 
nations of the two components were nearly identical, and I used 
the temperature as criterion instead. 

The inversions of the same spot on the two days agree that 
the more inclined component is weaker by around 0.5 kG in the 
innermost penumbra, but of same strength as the more vertical 



background component at the outer spot boundary due to a less 
steep decrease of fiel d strength with distan ce from the spot cen- 
ter, lik e also found bv lBorrero et al.l d2004l) or lBellot Rubio et ail 
(2004). I find hot upstreams in the mid penumbra, which appear 
in the more inclined component. The hot upflows show how- 
ever more vertical fields in the flow channel component than on 
neighboring pixels. The more inclined field component shows 
high flow velocities up to 5 kms -1 throughout the whole penum- 
bra, which I find to be roughly aligned with the field direction. 
At the outer spot boundary, the flow channels on average bend 
slightly downwards to return to the surface (or submerge below 
it), whereas the background component never exceeds an incli- 
nation of 60-70°. Both inverted sunspot maps show a change of 
the relative filling fraction of background component and flow 
channel component at around r/r spo , = 0.65, where the filling 
fraction of the flow channel component increases. 

To generate a geometrical model of the penumbral field 
topology, I integrated the surface inclination in radial direction. 
This approach is valid, if the field inclination (not field strength 
!) is independent of or only slightly dependent on depth in the 
solar atmosphere. This method can on the one hand be used as 
powerful visualization tool, as it allows to easily set geometry 
and all other field properties in context. On the other hand, inte- 
grating separately the inclination of background component and 
flow channel component, it allows to construct a 3-D model of 
the sunspot that is basically the best-fit atmosphere model for 
the observed profiles. In this 3-D model, the background com- 
ponent shows steeper ridges of enh anced field streng th, which 
maybe are identical to the "spines" o flLites et al.ldl993l) . The az- 
imuthally averaged flow channel component yields arched loops 
that cross the t = 1 surface at around r/r spo , = 0.65, where the 
filling fraction changes in favor of the flow channels. There is 
no obvious reason, why these two things happen at the same 
location: the inversion is done pixel-by-pixel without any infor- 
mation from neighboring pixels, whereas the integration uses all 
inclination values in radial direction. I thus conclude that it is 
highly probable that on average the flow channels do cross the 
surface at this radius. 

In general, the inversion results are i n agreement with the 
simula tions of the Moving Tube Model of ISchlichenmaier et al.l 
d 19981) in several aspects (geometry of background component 
and flow channels, radial variation of physical quantities). With 
a characteristic penumbral time scale of intensity variations be- 
low 30 min, and a depth of the penumbra of some Mm as sug- 
gested by the integrated inclination, I think that nothing found in 
the present observations would be in contradiction to a penum- 
bral heat transport by a series of hot, consecutively rising flow 
channels, as depicted in Fig. [15] An analysis of a time series of 
sunspot observations using a two-component inversion and the 
integration of the field inclination may be able to give a more di- 
rect proof if this scenario is actually happening in the penumbra. 
I will attempt this task as next step using either the data shown 
in Fig.[T2]or observations from later times t hat were taken with 
the he lp of adaptive optics system (e.g., ICabrera Solana et"aT1 
(12006k 
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Fig.A.l. A simple method to produce peculiar V profiles. The 
addition of the two profiles A and B in the top row yields the 
profiles in the bottom row. Left column: profile B equals -90 % 
of A and is located at the same wavelength. Right column: B 
equals -90 % of A, but now has been shifted in wavelength. 

Appendix A: Creation of multi-lobed profiles 

The Stokes V profiles in the neutral line on the limb side (Fig.[2]i 
show a complex multi-lobed structure. Although the profiles 
strongly differ from those originating from a simple atmosphere 
with constant magnetic field properties in the formation height 
of the spectral lines, only a few ingredients are necessary to re- 
produce their shapes to first order. Figure IA.1I shows a simple 
experiment: the regular profile from the center side (rightmost 
graph of Fig. [2]i is taken to be component A; a second com- 
ponent B is constructed as -90 % of A. Addition of these two 
profiles generates again a regular two-lobed profile with reduced 
amplitude (left column of Fig. [ATTT >. If component B is displaced 
in wavelength - corresponding to the Doppler shift induced by a 
flow field of 2 kms~' - the addition yields peculiar multi-lobed 
profiles like those observed in the neutral line of Stokes V. 

Appendix B: 2/3-D models 

B.I. 2-D: Integration of single column 

To follow the properties in a spatial cut without averaging, I 
took the inversion results along a single column of the 2-D maps 
(cf. Sect |4.4| i. The inclination of the cut can also be integrated in 
radial direction. In Fig. IB. Il l display the result in the same way 
as in Fig. [14] but this time with the temperature as color coding. 
The hot upstream is located at x = 9 Mm. The integrated curve 
was corrected to meet z=0 km at the outer sunspot boundary. 
Then I suggestively have shifted the integrated curve for x > 9 
Mm down by 150 km, to meet the r = 1 level at x = 9 Mm. As 
the formation height covers some hundred kilometers, the ge- 
ometry depicted still would be in agreement with the observed 
spectra, even if of course I have no proof that the second channel 
is below the first one. 





Fig. B.l. The integrated inclination values along the single col- 
umn marked in Fig. [8] Similar to Fig. [14] but using temperature 
for the color coding. The integrated curve has been shifted down 
at x = 7 Mm, where the hot upstream appears. 



B.2. 3-D models 



Figure IBT21 shows the full 3-D model for both observations from 
three different viewing directions: from above, the side, and 



slightly below the surface. The color code and display method 
correspond to that of Fig. QT] Two animations showing the 
sunspot model of 9.8.2003 are available in the online material. 
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Fig. B.2. The 3-D models from the integration of the inclinations to the surface averaged over 4 degree bins. {Left column): data 
of 9th of August. {Right column): data of 7th of August. {Top to bottom): view from above, side view, view from below. The color 
coding is identical to Fig.QT] The plot of the flow channels was switched off for the view from below. 



Appendix C: Analysis of DOT time series 

To quantify the amount of dynamics in the penumbra, I used a 2 
hours time series of speckle-reconstructed images in the G-band 
from the DOT telescope, taken on 9.8.2003 after the map used 
for the present investigation (30° heliocentric angle). To derive 



the boundaries of the penumbra, I used the temporal average of 
the time series (cf. Fig. IC.ll l. To avoid counting brightenings at 
the outer white-light boundary as penumbral grains (PGs), I fur- 
ther restricted the area to all points inside an ellipse centered on 
the spot, but with smaller radius than the average penumbra. I 
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Fig. C.l. Overview of the DOT data set. (Top row): single in- 
tensity image (left) with corresponding mask of brightenings 
(right). (Bottom row): Temporal averages over the time series. 
The black contours outline the area considered to be the penum- 
bra. Tick marks are in arcsec. 



created a mask of all points inside the penumbra above a thresh- 
old of 1 .2 times the continuum intensity outside the spot for each 
of the 214 images of the time series. The number and area of the 
PGs were taken from the statistics of the individual masks. 

On a single radial cut starting at the center and ending at the 
outer penumbral boundary, usually two or three individual PGs 
can be detected, even if they not always all exceed the threshold 
used. The temporal average of the masks clearly shows the de- 
pendence of the PGs on the azimuthal position inside the spot: on 
the limb side, the PGs appear roundish and have a much higher 
frequency than for azimuths of 90°, respectively, 270°, where 
they usually appear as thin streaks. On the center side, they are 
almost completely missing. As in a single image, usually two or 
more PGs can be seen for any radial cut in the temporal average. 
The intensity of all points inside the boundaries defined above 
during the time series yielded the histogram of Fig. IC.2I which 
shows that around 10 % of all points are on average brighter than 
the outside continuum intensity. 

All following statistics and quantities have been derived 
from the full time series in the area defined to belong to the 
penumbra. On average, in each image of the time series I found 
around 90 penumbral grains with an intensity above 1.2 times 
the continuum intensity outside the spot. The average size of the 
PGs was 16 pixels, which would correspond to a square area of 
200x200 km 2 . They covered an area fraction of around 3 % of 
the whole penumbra. Points with an intensity above a continuum 
intensity of unity covered 10 % of the full area. 

The intensity patterns in the penumbra disappear after a life- 
time of around 30 min; the autocorrelation drops to zero around 
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Fig. C.2. Average histogram of penumbral intensities in the time 
series. 
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Fig. C.3. Autocorrelation functions of the intensity in penumbra 
and a quiet Sun reference area. 

that time (cf. Fig. IC.3l ). For comparison, I used a granulation area 
outside the spot, which shows a characteristic correlation time 
of below 15 minutes. I used two types of autocorrelation with 
the same result: autocorrelation of the full images (restricted to 
the penumbra), and pixel-by-pixel autocorrelation, where only 
the intensity of one spatial position with time was used in the 
derivation of the autocorrelation. The number of PGs is however 
constant to a high degree. Thus, I conclude that the rate of the 
appearance of PGs should also have a characteristic timescale of 
around 30 min. 



